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Abstract 

Introduction: SNPs rs2981 582 and rs2981 578, located in a linkage disequilibrium block (LD block) within intron 2 
of the fibroblast growth factor receptor 2 gene {FGFR2), are associated with a mildly increased breast cancer risk. 
Allele-specific regulation of FGFR2 mRNA expression has been reported previously, but the molecular basis for the 
association of these variants with breast cancer has remained elusive to date. 

Methods: mRNA levels of FGFR2 and three fibroblast growth factor genes {FGFs) were measured in primary fibroblast 
and epithelial cell cultures from 98 breast cancer patients and correlated to their rs2981 578 genotype. The 
phosphorylation levels of downstream FGFR2 targets, FGF receptor substrate 2a (FRS2a) and extracellular signal- 
regulated kinases 1 and 2 (ERK1/2), were quantified in skin fibroblasts exposed to FGF2. Immunohistochemical markers 
for angiogenesis and lymphocytic infiltrate were semiquantitatively assessed in 25 breast tumors. 

Results: The risk allele of rs2981 578 was associated with increased FGFR2 mRNA levels in skin fibroblasts, but not in 
skin epithelial cell cultures. FGFR2 mRNA levels in skin fibroblasts and breast fibroblasts correlated strongly in the 
patients from whom both cultures were available. Tumor-derived fibroblasts expressed, on average, eight times 
more FGFR2 mRNA than the corresponding fibroblasts from normal breast tissue. Fibroblasts with higher FGFR2 
mRNA expression showed more FRS2a and ERK1/2 phosphorylation after exposure to FGF2. In fibroblasts, higher 
FGFR2 expression correlated with higher FGF10 expression. In 25 breast tumors, no associations between breast 
tumor characteristics and fibroblast FGFR2 mRNA levels were found. 

Conclusions: The influence of rs2981 578 genotypes on FGFR2 mRNA expression levels is cell type-dependent. 
Expression differences correlated well with signaling levels of the FGFR2 pathway. Our results suggest that the 
increased breast cancer risk associated with SNP rs2981 578 is due to increased FGFR2 signaling activity in stromal 
fibroblasts, possibly also involving paracrine FGF10 signaling. 



Introduction 

Several genome-wide association studies have shown 
that the minor allele of SNP rs2981582, located in 
intron 2 of the fibroblast growth factor receptor 2 gene 
(FGFR2), is associated with increased breast cancer risk 
[1-6]. The odds ratios (ORs) for this SNP are 1.23 in 
heterozygotes and 1.63 in homozygotes for the minor 
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allele [1,2], but attributable risk is high because of the 
high frequency of the risk allele population. Fine-scale 
genetic mapping and resequencing of the region sur- 
rounding rs2981582 resulted in the identification of up 
to eight variants in a linkage disequilibrium block (LD 
block) within intron 2 of FGFR2 most strongly asso- 
ciated with increased breast cancer risk, including SNP 
rs2981578 [1,5,7]. The location of this LD block suggests 
that these variations somehow modify the functioning of 
FGFR2. 
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Meyer et al. [7] have shown that two SNPs within this 
LD block, one of them being rs2981578, alter the DNA 
binding affinity of octamer-binding transcription factor 
1 (Oct-1), runt-related transcription factor 2 (Runx2) 
and CCAAT/enhancer binding protein P (C/EBPp 1 ). 
Accordingly, increased expression of FGFR2 mRNA was 
observed in total RNA isolated from breast tumors of 
patients homozygous for the risk allele as compared to 
homozygotes for the major allele [7]. Paradoxically, Sun 
et al. [8] recently reported decreased expression of 
FGFR2 mRNA in normal breast tissue of homozygotes 
for the risk allele. 

FGFR2 is one of five fibroblast growth factor (FGF) recep- 
tors known in humans to be involved in various signaling 
pathways that regulate processes such as cell growth, apop- 
tosis and differentiation. Two isoforms, FGFR2-IIIb and 
FGFR2-IIIc, are the result of mutually exclusive alternative 
splicing of exon 9 or 10 of FGFR2. Isoform Illb is present 
on epithelial cells and binds ligands FGF3, FGF7, FGF10 
and FGF22 and isoform IIIc is present on mesenchymal 
cells and binds FGF2, FGF4, FGF6, FGF9, FGF17 and 
FGF18 [9]. Binding of a ligand to the receptor can activate 
several signaling pathways, including the mitogen-activated 
protein kinase (MAPK) pathway [10,11]. 

Downregulation of FGFR2 protein has been reported in 
up to 67% of breast tumors [12], whereas amplification of 
FGFR2 and upregulation of FGFR2 mRNA expression 
have been reported in less than 10% of breast tumors 
[10,13]. Somatic FGFR2 mutations are rare in breast can- 
cer [14]. A switch from the Mb to the IIIc isoform in 
tumor cells, resulting in activation of the receptor by 
other FGFs, has been reported in a subset of prostate 
cancers and in a few breast cancer cell lines [15-17]. 

To further unravel the mechanisms by which the SNPs 
in intron 2 of FGFR2 increase breast cancer risk, and to 
address the heterogeneous cellular composition of breast 
tumors, we studied the expression of FGFR2 mRNA in 
relation to the rs2981578 genotype in fibroblasts and 
epithelial cells cultured from breast tissue. We also com- 
pared the FGFR2 mRNA expression in fibroblasts derived 
from normal skin tissue, normal breast tissue and breast 
tumor tissue. Furthermore, we explored the functional 
implications of different levels of FGFR2 expression at the 
cellular level by studying the downstream MAPK pathway. 
Finally, we compared histological characteristics between 
tumors of patients with high and low FGFR2 levels. 

Materials and methods 

Patients 

Breast tissue was obtained from 98 women who had 
undergone breast cancer-related surgery at the Leiden 
University Medical Center (LUMC) from February 2006 
until December 2010 (patient demographics are given in 
Table 1). According to the Dutch Medical Treatment 



Act, anonymized samples obtained during medical treat- 
ment may be used in medical research if the patient 
does not object to this "secondary use" (opt-out system). 
Samples were coded in such a way that identifying infor- 
mation was not accessible by the researcher, but sam- 
ples still could be linked to pathology reports by 
authorized clinicians. The use of residual material for 
this study was approved by the Department of Pathology 
of the LUMC (code OP 25-06). The study protocol was 
compliant with the Code of Conduct issued by the 
Dutch Federation of Medical Scientific Societies. 

Tissue culture 

Fibroblasts and epithelial cells were cultured separately 
from skin tissue, normal breast tissue or tumor tissue by 
cutting the tissue patches into 2-mm x 2-mm pieces. 
Tissue patches were allowed to attach to the culture 
flask for 30 minutes, then we added DMEM/F12 med- 
ium (1:1 dilution; GIBCO/Invitrogen, Breda, the Nether- 
lands) supplemented with 20% FCS (Bodinco, Alkmaar, 
The Netherlands) and antibiotics. Cells were grown at 
37°C in 5% C0 2 . Each fibroblast culture was expanded 
in medium supplemented with 10% FCS and antibiotics. 
Epithelial cells were grown in HuMEC Ready Medium 
(GIBCO/Invitrogen) or in DermaLife K Serum-Free Ker- 
atinocyte Culture Medium (Lifeline Cell Technology, 
Walkersville, MD, USA) (see also Table 2). The mor- 
phology of primary cultures was in agreement with the 
morphology of established epithelial and mesenchymal 
cell lines (Figure SI in Additional file 1). 

Genotyping 

Genomic DNA was isolated using the Wizard Genomic 
DNA Purification Kit (Promega, Leiden, The Netherlands) 
according to the manufacturer's instructions. Four SNPs 
in intron 2 of FGFR2 were genotyped (Table 3). 

Table 1 Characteristics of the patient cohort 3 

Characteristics Data 

Total number of patients 98 

BRCA1/2 mutation carriers, n (%) 10 (10%) 
Diagnosis 

Mean age at diagnosis' 5 , years 55 (±14) 
(±SD) 

Invasive breast cancer, n 83 

In situ breast cancer, n 8 

Prophylactic surgery, no cancer, 7 
n 

Tumor characteristics (n = 91) 

ER status 54 positive, 20 negative, 17 

unknown 

a SD: standard deviation; ER: estrogen receptor. b Age at diagnosis of breast 
cancer. The seven patients who underwent prophylactic surgery were 
excluded from this analysis. 
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Table 2 Culture medium used to culture skin epithelial 



cells 3 


Culture ID 


Medium used 


46 


HuMEC 


47 


HuMEC 


52 


HuMEC 


54 


HuMEC 


55 


HuMEC 


57 


HuMEC 


58 


HuMEC 


59 


HuMEC 


64 


HuMEC 


65 


HuMEC 


69 


DermaLife 


70 


DermaLife 


77 


DermaLife 


83 


HuMEC 


100 


HuMEC 


1 03 


L/tri I I id i_i I c 


132 


DermaLife 


140 


DermaLife 


143 


DermaLife 


152 


DermaLife 


157 


DermaLife 


163 


DermaLife 


170 


DermaLife 


189 


DermaLife 


209 


DermaLife 



a Culture ID refers to the unique number assigned to a patient during the 
process of sample anonymization; HuMEC = HuMEC Ready Medium (GIBCO/ 
Invitrogen, Breda, The Netherlands}; DermaLife = DermaLife K Serum-Free 
Keratinocyte Culture Medium (Lifeline Cell Technology, Walkersville, MD, USA). 



Quantitative real-time PCR 

Total RNA was isolated from exponentially growing 
fibroblasts or epithelial cells using NucleoSpin RNA II 
Kit (Clontech, St-Germain-en-Laye, France) according 



to the manufacturer's instructions. During isolation, a 
30-minute DNase treatment was performed to ensure 
complete degradation of contaminating DNA. 

A detailed description of the RT reaction and quantita- 
tive PCR (qPCR) conditions can be found in the supple- 
mentary methods section (Additional file 2). Briefly, qPCR 
was performed to analyze the expression of a specific gene 
(Table 4) using the Bio-Rad iCycler thermal cycler and the 
iQ SYBR Green Supermix (Bio-Rad Laboratories, Her- 
cules, CA, USA) according to the manufacturer's recom- 
mendations. Each gene was measured at least twice in 
each sample. As reference genes, the TATA box binding 
protein (TBP) and heterogeneous nuclear ribonucleopro- 
tein M (HNRPM or HNRNPM) (fibroblasts) genes or sig- 
nal recognition particle receptor {SRPR), HNRPM and 
TBP (epithelial cells) genes were used. 

Fragment analysis 

Isoforms FGFR2-IIIb and FGFR2-IIIc differ by the alter- 
native splicing in of either exon 9 or exon 10 into the 
mature mRNA, resulting in a difference of 3 bp in 
length of the mature mRNA. A PCR fragment of 297 bp 
for FGFR2-Ulb or 300 bp for FGFR2-lllc mRNA was 
amplified by PCR with 38 cycles and an annealing tem- 
perature of 70°C using the 6-carboxyfluorescein-labeled 
forward primer 5'-GTGGAAAAGAACGGCAGTAAA- 
TACG-3' and reverse primer 5'-CACCATACAGGC- 
GATTAAGAAGACC-3' located in exons 8 and 11. 

Fragment analysis was performed on cDNA from 14 
skin epithelial cell cultures and 52 skin fibroblast cul- 
tures using the 3730x1 DNA Analyzer (Applied Biosys- 
tems, Carlsbad, CA, USA) and the GeneScan 400HD 
ROX Size Standard (Applied Biosystems). 

Measuring FGFR2 protein levels 

FGFR2 protein levels were analyzed in seven exponen- 
tially growing fibroblast cultures with high (n = 3) or 



Table 3 Primers used to genotype four SNPs in FGFR2* 



SNP 


Primer 


Primer sequence (5' to 3') 


Product 
length 


Annealing 
temperature 


Further analysis 




rs 10736303 


F 
R 


AGGGACAAATACTCCGCACA 
AGCCATCCAGCATGmCTC 


405 bp 


50°C 


Sanger sequencing 




rs2981578 


F 
R 


TGACTCTTCAAAG 1 1 IGI 1 IGI 1 1 1 

GAGGAAAGGTTCCCCACACT 


295 bp 


50°C 


Restriction enzyme Aa\ (New England 
MA, USA) 


3iolabs, Ipswich, 


rs2981582 


F 
R 


AGCTCAGCTFACCCCAGACA 
CGTGAGCCAAGCCTCTACTT 


215 bp 


58°C 


Restriction enzyme Acil (New England 
MA, USA) 


3iolabs, Ipswich, 


rs7895676 


F 
R 


CAGGTGCGGTGGCTCATGTC 
GACTTCAATGGCGGGACTCC 


345 bp 


67°C 


Sanger sequencing 





a FGFR2: fibroblast growth factor receptor 2 gene; F: forward primer; R: reverse primer; bp: base pair. For each SNP, the forward and reverse primer sequences, the 
length of the resulting PCR product and the annealing temperature used during PCR are shown. Also shown is the technique that was used on the resulting PCR 
fragment to genotype the SNPs. 
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Table 4 Primers used during real-time PCR analyses 3 


Gene and 
(GenelD) b 


Sequence 
ID 


Primer 


Primer sequence (5' to 3') 


Product 
length 


Distance from poly(A) 
tail c 


Exon-exon 
boundary 


FGF2 


NM 002006.4 


F 


ACG I GGAGAG I GG I I I I I GCGGA 


91 bp 


1,731 bp 


No 


(2247) 




R 


GGTGCCACGTGAGAGCAGAGC 








FGF7 


NM 002009.3 


F 


CTCAACGGCAAGTTTCGCTCCG I I I I C 


80 bp 


1,701 bp 


No 


(2252) 




R 


GCCTTCCAGGATTTGCTGGCCC 








FGF10 


NM 004465.1 


F 


TCTTCTTCCTCCTGCTTCTCCTCTCC 


148 bp 


325 bp 


No 


(2255) 




R 


TCCCGCTGACCTTCCCG^CTTCTC 








FGFR2 d 


NM 000141.4 


F 


GTCAGTGAGAACAGTAACAACAAG 


192 bp 


3,368 bp 


Yes 


(2263) 




R 


GTAGCCTCCAATGCGATGC 








GADPH 


NM 002046.3 


F 


TTCCAGGAGCGAGATCCCT 


175 bp 


805 bp 


Yes 


(2597) 




R 


CACCCATGACGAACATGGG 








HMBS 


NM 000190.3 


F 


CTGGTAACGGCAATGCGGCT 


338 bp 


1010 bp 


Yes 


(3145) 




R 


GCAGATGGCTCCGATGGTGA 








HNRPM 


NM_005968.3 


F 


GAGGCCATGCTCCTGGG 


85 bp 


440 bp 


Yes 


(4670) 




R 


TTTAGCATCTTCCATGTGAAATGG 








HPRF1 


NM_000 194.2 


F 


TGACACTGGCAAAACAATGCA 


94 bp 


742 bp 


Yes 


(3251) 




R 


GGTCCTTTTGACCAGCAAGCT 








SRPR 


NM_003139.2 


F 


CATTGCTTFTGCACGTAACCAA 


70 bp 


1,308 bp 


Yes 


(6734) 




R 


ATTGTCTTGCATGGGGCG 








FBP 


NMJXB 194.3 


F 


CACGAACCACGGCACTGATT 


89 bp 


905 bp 


Yes 


(6908) 




R 


mTCTTGCTGGCAGTCTGGAC 









a FGF: fibroblast growth factor gene; FGFR: fibroblast growth factor receptor gene; GAPDH: glyceraldehyde 3-phosphate dehydrogenase gene; HMBS: 
hydroxymethylbilane synthase gene; HNRPM: heterogeneous nuclear ribonucleoprotein M gene; HPRT1: hypoxanthine phosphoribosyltransferase 1 gene; SRPR: 
signal recognition particle receptor gene; TBP: TATA box binding protein gene; F: forward primer; R: reverse primer; bp: base pair. For each gene, the forward and 
reverse primer sequences, the length of the resulting PCR product (using cDNA as template) and the distance from the poly(A) tail are shown, as is the 
information about whether the PCR product includes an exon-exon boundary. b GenelD from EntrezGene accessed 2 June 2010. c Distance from 3' end of the PCR 
product to 5' beginning of the poly(A) tail. d Primer set for FGFR2 is in a region of FGFR2 that is present in both FGFR2 1Mb and FGFR2 lllc. 



low (n = 4) FGFR2 mRNA levels. Total protein isolation 
and Western blot analysis are described in detail in the 
supplementary methods section (Additional file 2). Each 
sample was measured twice on separate gels. The pri- 
mary and secondary antibodies used are given in Table 
5. The Odyssey Infrared Imaging System and Odyssey 
3.0 software (LI-COR Biotechnology, Cambridge, UK) 
were used to visualize and identify the bands. 

Analyzing the phosphorylation of downstream FGFR2 
targets 

Five skin fibroblast cultures with high FGFR2 mRNA 
expression and five skin fibroblast cultures with low 
FGFR2 mRNA expression were randomly selected. In 
each experiment, one culture with high mRNA expres- 
sion and one culture with low mRNA expression were 
compared. In each well of a six-well plate, 5 x 10 fibro- 
blasts from the full cell culture flasks were seeded and 
10% FCS medium was added. The cells were allowed to 
attach to the well bottom for 24 hours and were then 
starved in 0.5% FCS medium for 12 hours. Cells were 
subsequently exposed to 0 or 20 ng of FGF2 (basic 
fibroblast growth factor, 9952; Cell Signaling Technol- 
ogy, Danvers, MA, USA) in 0.5 mL of 0.5% FCS medium 
for 10 minutes. This length of time was chosen because 



in an initial experiment phosphorylation levels peaked at 
about 10 minutes (Figure S2 in Additional file 3). 

Total protein was then isolated, and Western blot ana- 
lysis was performed to assess the phosphorylation of 
downstream targets of FGFR2 as described in more 
detail in the supplementary methods section (Additional 
file 2). The primary and secondary antibodies that were 
used are given in Table 5. Phosphorylation levels after 
10 minutes of exposure were normalized to phosphory- 
lation levels before exposure. 

Determining the stroma percentage within the primary 
tumor 

Routine H & E-stained histological tumor sections were 
obtained from the Department of Pathology for 50 of 
the 68 patients from whom skin fibroblast cultures were 
available. The stroma percentage within the primary 
tumor was estimated by two researchers as previously 
described [18]. 

Immunohistochemistry 

For the immunohistochemical analyses, patients were 
selected for being either homozygous for the risk allele 
and having high FGFR2 mRNA expression (n = 11) or 
for being homozygous for the major allele and having 
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Table 5 Antibodies used during Western blot analyses 3 



Antigen targeted 



Antibody (catalog number, manufacturer) 



Clonality 



Dilution 



Origin 



Primary antibodies 
FGFR2 
a-tubulin 
ERK1/2 

Phosphorylated ERK1/2 
FRS2a 

Phosphorylated FRS2a 
Secondary antibodies 
Mouse IgG 
Rabbit IgG 



Bek C-17 (sc-122, SZ) 
Anti-a-tubulin (T6199, SA) 
p44/42 MAPK antibody (9102, CS) 
Phospho-p44/42 MAPK antibody (4377, CS) 
FRS2 H-91 (sc-8318, SZ) 
Phospho-FRS2-a antibody (3864, CS) 

IRDye 680 anti-mouse IgG (926-32220, LC) 
IRDye 800CW anti-rabbit IgG (926-3221 1, LC) 



Polyclona 
Monoclonal 

Polyclonal 
Monoclonal 

Polyclonal 

Polyclona 

Polyclonal 
Polyclonal 



in 500 
in 8,000 
in 100 
in 200 
1 in 50 
1 in 50 

in 7,500 
in 5,000 



Rabbit 
Mouse 
Rabbit 
Rabbit 
Rabbit 
Rabbit 

Goat 
Goat 



a FGFR2: fibroblast growth factor receptor 2 gene; ERK1/2: extracellular signal-regulated kinases 1 and 2; FRS2a: fibroblast growth factor receptor substrate 2a; 
IgG: immunoglobulin G; MAPK: mitogen-activated protein kinase. Manufacturer names are abbreviated as follows: SZ: Santa Cruz Biotechnology, Santa Cruz, CA, 
USA; SA: Sigma-Aldrich, Zwijndrecht, The Netherlands; CS: Cell Signaling Technology, Danvers, MA, USA; LC, LI-COR Biosciences, Cambridge, UK. For each 
antibody, the target, manufacturer, clonality, dilution used during Western blot analysis and origin are given. 



low FGFR2 mRNA expression (n = 14). Tissue sections 
were stained using antibodies to detect the presence of 
microvessels (CD31) and T lymphocytes (CD3). Micro- 
vascular density was measured by two independent 
observers who counted all CD31 -stained blood vessels 
in 10 high-power fields (HPFs) randomly selected within 
the tumor region. The average number of blood vessels 
per HPF was then calculated and used in subsequent 
analysis. For CD3, slides were scanned using the Pan- 
noramic MIDI digital scanner (3DHISTECH Ltd., Buda- 
pest, Hungary) and then viewed at xlO magnification 
using the MIRAX Viewer 1.11 (3DHISTECH Ltd.). For 
each tumor, four to six representative areas of 1,800 um 
x 1,200 um were selected, and the tumor and surround- 
ing stroma were separately analyzed using ImageJ soft- 
ware and the color deconvolution plug-in [19]. The area 
covered by the antibody was measured and normalized 
by dividing it by the area covered by nuclei. 

Statistical analysis 

The distribution of the SNP genotypes was studied 
using an Excel-based Hardy-Weinberg equilibrium 
(HWE) calculator [20]. For all other statistical analyses, 
we used PASW Statistics version 17.0 software (SPSS 
Inc., Nieuwegein, The Netherlands). 

Results 

A total of 68 skin fibroblast cultures, 44 breast fibroblast 
cultures, 11 tumor-derived fibroblast cultures and 25 
skin epithelial cell cultures were established from 98 
patients who had undergone breast cancer-related sur- 
gery at LUMC. We genotyped four SNPs in intron 2 of 
FGFR2 (Table 6). All SNPs were in HWE. SNP 
rs2981578 was chosen to represent all four SNPs 
because of its reported effect on the binding of tran- 
scription factor Runx2 and because it showed strong LD 
with the other three SNPs including rs2981582. 



FGFR2 mRNA expression in fibroblasts and epithelial cells 

A significant correlation between the rs2981578 geno- 
type and FGFR2 mRNA expression level as measured by 
qRT-PCR was present in the 68 skin fibroblast cultures 
(Figure 1A) (P = 0.02; one-way analysis of variance 
(ANOVA)). Although expression levels varied widely 
within each genotype group, a dosage-dependent effect 
was seen, with the highest average expression found 
among homozygotes for the minor allele (the risk allele), 
followed by heterozygotes, and the lowest average 
expression found among homozygotes for the major 
allele. In contrast, among 25 skin epithelial cell cultures, 
there were no statistically significant differences in 
FGFR2 mRNA expression levels between the three gen- 
otype groups, although the average levels of FGFR2 
mRNA were slightly lower in carriers of one or two 
copies of the risk allele (Figure IB) (P = 0.732; one-way 
ANOVA). 



Table 6 Genotype distribution of four SNPs in FGFR2" 



rs2981582 rs2981578 rsl 0736303 rs7895676 Number of 

patients 



CT 


AG 


AG 


CT 


29 


TT 


GG 


GG 


CC 


22 


CC 


AA 


AA 


TT 


21 


CT 


GG 


GG 


CC 


7 


CC 


AG 


AG 


CT 


5 


TT 


AA 


AA 


TT 


1 


Total 








85 b 



a FGFR2: fibroblast growth factor receptor 2 gene. The number of patients who 
carry each combination of the four genotypes is shown in the right-hand 
column. Although we have no haplotype information, the genotype 
distribution shows the strong linkage disequilibrium present within this 
region. The alleles associated with increased breast cancer risk are rs2981582: 
T; rs2981578: G; rs10736303: G; and rs7895676: C. b For the remaining13 
patients, only rs2981578 was genotyped; of these patients, 5 carried AA, 5 
carried AG and 3 carried GG. 
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minor 
homozygotes 



rs2981578 genotype 



1 .5—1 



S I- 



D) 

§0,5- 



«> o- 

< 

z 

DC 

E-0,5" 



O -1- 



-1,5- 



minor 
homozygotes 



1 

hetero zygotes 



major 
homozygotes 



rs2981578 genotype 



Figure 1 The relationship between FGFR2 mRNA expression 
and the rs2981578 genotype in skin fibroblast and skin 
epithelial cell cultures. (A) A significant correlation between the 
rs2981 578 genotype and FGFR2 mRNA expression level as measured 
by quantitative real-time PCR was present in the 68 skin fibroblast 
cultures (P = 0.02; one-way ANOVA). The expression levels were 
log 2 -transformed and normalized to HNRPM and TBP expression. (B) 
FGFR2 mRNA expression and rs2981 578 genotypes in 25 skin 
epithelial cell cultures (P = 0.73; one-way ANOVA). The expression 
levels were log 2 -transformed and normalized to HNRPM, SRPR and 
TBP expression. Each dot represents the expression level in one 
patient. 



FGFR2 mRNA expression in fibroblasts derived from 
different tissues 

To address whether FGFR2 expression levels in skin 
fibroblasts correlated to those in breast fibroblasts, we 
cultured fibroblasts from the skin tissue as well as the 
normal breast tissue from 22 breast cancer patients. 
There was a significant correlation between the FGFR2 
mRNA levels in the skin fibroblasts and the breast 



fibroblasts (Figure 2A) (Pearson's correlation coefficient 
= 0.64, P = 0.001). Expression levels were on average 2.8 
times higher in the breast fibroblasts than in the skin 
fibroblasts {P = 0.001; paired f-test). 

For 11 of these 22 patients, we also successfully grew 
fibroblasts from the tumor tissue. FGFR2 mRNA expres- 
sion in tumor-derived fibroblasts was on average eight 
times higher than in breast fibroblasts and sixteen times 
higher than in skin fibroblasts (Figure 2B) (P = 3 x 10~ 4 
and P = 2 x 10' , respectively; paired £-tests). 




-6-4-2 0 2 

FGFR2 mRNA level (2log scale) in skin fibroblasts 




skin normal breast tumor 

Origin of the cultured fibroblasts 

Figure 2 FGFR2 mRNA expression levels in fibroblasts cultured 
from different locations. (A) We found a significant correlation 
between FGFR2 mRNA expression in skin fibroblasts and normal 
breast fibroblasts from the same patient (n = 22). The FGFR2 
expression levels were log 2 -transformed and normalized to HNRPM 
and TBP expression. (B) FGFR2 mRNA expression in tumor-derived 
fibroblasts is significantly higher than in fibroblasts cultured from 
normal breast tissue. The FGFR2 expression levels were log 2 - 
transformed and normalized to HNRPM and TBP expression. Each 
dot represents FGFR2 expression in one patient. The lines connect 
the fibroblast cultures from one patient. The horizontal bars shows 
the mean expression for all fibroblast cultures from that origin (skin 
-3.8 (±1.9 SD), normal breast -2.9 (±1.4 SD), tumor tissue 0.2 (±1.5 
SD)). *P < 0.001. 
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FGFR2 mRNA isoform expression analysis in fibroblasts 
and epithelial cells 

Since an isoform switch from FGFR2 Illb to FGFR2 IIIc 
has been reported in a few breast cancer cell lines, and 
since this switch may contribute to tumor development, 
we explored whether rs2981578 genotype influences the 
alternative splicing of isoforms of FGFR2. The expected 
FGFR2-Yllb was detected in 14 epithelial cell cultures, 
and the expected FGFR2-l\\c was detected in 48 of 52 
fibroblast cultures. In the remaining four fibroblast sam- 
ples, FGf7?2-IIIb was present at a level approximately 
10% of that of FGfl?2-IIIc, but the presence or absence 
of FGi\R2-IIIb was not associated with rs2981578 geno- 
type (P = 0.32; Pearson's % 2 test). 

Correlation between FGFR2 and FGF10 mRNA expression 

To explore a possible paracrine relationship between 
fibroblasts and epithelial cells, the mRNA expression 
levels of FGF2, FGF7 and FGF10 were measured in the 
68 skin fibroblast cultures and compared to their FGFR2 
mRNA expression. FGF2 mRNA expression was not 
correlated with FGFR2 expression (Figure 3C) (Pearson's 
correlation coefficient = 0.22, P = 0.07). The correlation 
between FGFR2 and FGF7 mRNA levels in skin fibro- 
blasts did not remain significant after correction for 
multiple testing (Figure 3B) (Pearson's correlation coeffi- 
cient = 0.28, P = 0.02). 

FGFR2 and FGF10 expression levels were loosely but 
significantly correlated in the skin fibroblast cultures 
(Figure 3A) (Pearson's correlation coefficient = 0.60, P = 
8 x 10" 8 ). There was also a suggestive correlation 
between FGF10 expression and rs2981578 genotype in 
these cultures (P = 0.06; one-way ANOVA) (Figure S3 
in Additional file 4). In 44 fibroblast cultures from nor- 
mal breast tissue, there was a suggestive correlation 
between FGFR2 and FGF10 mRNA levels (Spearman's p 
= 0.25, P = 0.11) (Figure S4 in Additional file 5). 

FGFR2 protein expression 

In seven fibroblast cultures, FGFR2 protein levels were 
quantified by Western blot analysis and compared with 
the FGFR2 mRNA levels. The FGFR2 antibody that was 
used detected a full-length FGFR2 recombinant protein, 
but absolute protein expression levels in fibroblasts were 
very low, giving poor signal-to-noise ratios. Possibly as a 
result of this, we were unable to demonstrate a correla- 
tion between FGFR2 mRNA and protein levels (Spear- 
man's p = -0.14, P = 0.76) (Figure S5 in Additional file 6). 

Phosphorylation of downstream FGFR2 targets 

To investigate whether the differences in FGFR2 mRNA 
expression levels correspond with different levels of 
FGFR2 pathway activity, we compared the phosphoryla- 
tion of downstream targets FGF receptor substrate 2a 
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Figure 3 Correlation of expression levels of FGFR2 and FGF10 
mRNA in 68 skin fibroblast cultures No correlation was found for 
FGFR2 and FGF2 mRNA or FGFR2 and FGF7 mRNA levels. The data 
were normalized to HNRPM and TBP and log 2 -transformed. Each dot 
represents the expression levels from one patient. (A) FGFR2 and 
FGF10 (P = 8 x 10* Pearson's correlation). (B) FGFR2 and FGF7 (P = 
0.02; Pearson's correlation). (C) FGFR2 and FGF2 (P = 0.07; Pearson's 
correlation). 
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(FRS2a) and extracellular signal-regulated kinases 1 and 
2 (ERK1/2) after 10 minutes of exposure to ligand FGF2 
in skin fibroblast cultures with high and low FGFR2 
mRNA levels. Duplicate measurements correlated well 
in all experiments (Pearson's correlation coefficient = 
0.92, P = 8 x 10" 65 ). 

Phosphorylation of FRS2a was measured in four fibro- 
blast cultures, two with high and two with low FGFR2 
expression. The sample with higher FGFR2 expression 
showed a larger increase in phosphorylation of FRS2a 
after FGF2 was added compared to the sample with 
lower FGFR2 expression (Figure 4) (P = 0.18; Wilcoxon 
signed-rank test). 

Phosphorylation of ERK1/2 increased in all skin fibro- 
blast cultures after exposure to FGF2, and this effect 
was again more pronounced in the sample with higher 
FGFR2 expression in all five experiments (Figure 4) (P = 
0.04; Wilcoxon signed-rank test). 

FGFR2 mRNA level and histological breast tumor 
characteristics 

We explored tumor parameters in patients with differ- 
ent FGFR2 mRNA levels in their fibroblasts. Since the 
percentage of stroma in the primary tumor may vary 
considerably [18] and FGFR2 signaling in the stroma 
could influence cell proliferation, we explored whether 
patients with higher FGFR2 mRNA expression in fibro- 
blasts had a higher stroma percentage. In the 50 patients 



analyzed, no correlation was present between FGFR2 
mRNA levels in their fibroblasts and the percentage of 
stroma within the tumor (Spearman's p = 0.17, P = 
0.23). Unfortunately, not enough tumor tissue was avail- 
able from patients from whom we had obtained skin 
epithelial cells. 

Because FGF signaling in cancer has been associated 
with induction of angiogenesis and lymphocytic infiltra- 
tion [21], the number of microvessels was determined 
and T-cell infiltrate was quantified in breast tumors 
from patients whose fibroblasts had relatively high (n = 
11) or relatively low (n = 14) FGFR2 expression. Inter- 
observer agreement regarding the number of blood ves- 
sels per HPF was high (Pearson's correlation coefficient 
= 0.89, P = 1 x 10 s ), but no correlation with FGFR2 
mRNA levels was present (Pearson's correlation coeffi- 
cient = 0.04, P = 0.86). No correlations between the 
amount of T-cell infiltrate and FGFR2 mRNA levels in 
the fibroblasts were detected either (T cells in stroma: 
Pearson's correlation coefficient = 0.04, P = 0.86; T cells 
in tumor: Pearson's correlation coefficient = 0.15, P = 
0.52). 

Finally, since the SNPs in intron 2 of FGFR2 are 
mainly associated with estrogen receptor-positive (ER- 
positive) breast tumors [1,22], we studied a possible cor- 
relation between FGFR2 mRNA levels and ER status of 
the tumor in 50 patients, but could detect none (P = 
0.60; logistic regression). 




Figure 4 Response of skin fibroblast cultures with high and 
low FGFR2 mRNA levels to stimulation with ligand FGF2 

Phosphorylation of downstream targets FRS2a and ERK1/2 were 
measured before and after 10 minutes of stimulation with FGF2. 
The average difference of the increase in phosphorylation is shown, 
which was calculated by dividing the increase in samples with high 
FGFR2 mRNA levels by the increase in samples with low FGFR2 
mRNA levels. The 95% confidence intervals are also shown. 
Phosphorylation of FRS2a was measured in four fibroblast cultures: 
two with high and two with low FGFR2 mRNA levels. 
Phosphorylation of ERK1/2 was measured in 10 fibroblast cultures, 
five with high and five with low FGFR2 mRNA levels. 



Discussion 

Aberrant FGF signaling has been implicated in the 
pathogenesis of multiple types of cancer, including 
breast cancer [11,23]. Germline genetic variation in an 
LD region within intron 2 of FGFR2 has been associated 
with a modestly increased breast cancer risk. Fine-scale 
genetic mapping and functional analyses identified up to 
eight SNPs, including SNP rs2981578, as likely causal 
variants [1,5,7]. It was demonstrated that these SNPs 
influence mRNA expression levels of the FGFR2 gene. 
Despite this progress, substantial uncertainty remains 
regarding how subtle modulation of FGFR2 expression 
levels may modify breast cancer risk. In addition, oppo- 
site effects of intron 2 SNP genotypes on FGFR2 mRNA 
expression have been reported in tumor tissue and nor- 
mal breast tissue [7,8]. 

We found, on average, higher FGFR2 mRNA levels in 
skin fibroblast cultures of heterozygotes and homozy- 
gotes for the risk allele of SNP rs2981578 relative to 
homozygotes for the major allele. This is in agreement 
with measurements of FGFR2 mRNA expression in total 
RNA from breast tumor tissue [7]. However, no effect 
of genotype was observed in our skin epithelial cell cul- 
tures. Also, epithelial breast tumor cells express lower 
levels of FGFR2 than surrounding normal breast 
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epithelium [12]. Soluble factors secreted or recruited by 
tumor cells may induce expression changes in neighbor- 
ing normal cells, leading to gene expression differences 
between tumor-derived stromal fibroblasts and normal 
breast stroma [24]. In the 11 patients from whom we 
obtained fibroblasts from skin, normal breast tissue and 
breast tumor tissue, FGFR2 expression levels in the 
tumor-derived fibroblasts were consistently higher than 
in the fibroblasts from normal breast tissue. A similar 
finding was made recently in cell cultures from patients 
with esophageal cancer [25]. Unfortunately, our sample 
set was too small to conclusively establish whether 
expression levels of FGFR2 mRNA in tumor-derived 
fibroblasts are similarly associated with SNP genotypes 
in intron 2 as they are in skin-derived fibroblasts. Taken 
together, however, these data suggest that the associa- 
tion between genotype and FGFR2 expression observed 
by Meyer et al. [7] in total breast tumor homogenates 
was derived from cancer-associated fibroblasts rather 
than from the tumor (or epithelial) component. Because 
the transcription factors Oct-l/Runx2 and C/EBPp dis- 
play differential binding to these SNP alleles in epithelial 
breast tumor cells, it will be interesting to study this in 
other cell types constituting breast tumor stroma. 

Our results in fibroblast and epithelial cell cultures are 
at odds with those obtained by Sun et al. [8], who reported 
lower FGFR2 mRNA expression in normal breast tissue 
from homozygotes for the risk allele. However, the differ- 
ences reported by Sun et al. were small and the correlation 
was weak. In our epithelial cell cultures, there was a trend 
toward lower FGFR2 mRNA expression in carriers of one 
or two copies of the risk allele, but the differences were 
very small. The insufficient statistical power of our studies 
or the differences in experimental design (whole tissue 
analysis versus analysis of cultured epithelial cells) may 
underlie these different outcomes. 

The variation in FGFR2 mRNA expression between 
individuals with the same intron 2 genotype is wide in 
both fibroblasts and epithelial cells. Apparently, the cau- 
sal variants in intron 2 of FGFR2 only partly determine 
the FGFR2 mRNA expression level. We have shown that 
in skin fibroblasts, higher FGFR2 mRNA levels corre- 
spond to higher activity of the FGFR2 pathway upon sti- 
mulation by ligand FGF2, indicating that the allelic 
status of the intron 2 SNPs is functional at the FGF sig- 
naling level. Cells with higher FGFR2 expression prob- 
ably respond differently to the same fibroblast growth 
factor concentration in their microenvironment [26]. 
Thus, the overall activity of the FGFR2 pathway might 
be a more accurate predictor of breast cancer risk than 
rs2981578 or rs2981582 genotype status, although it is 
presently unclear through which cell type and in which 
developmental phase of the mammary gland this activity 
contributes most to this risk. 



Fibroblasts have a well-recognized role in the carcino- 
genic process as remodelers of the extracellular matrix 
in tumor stroma and as a source of paracrine growth 
factors that influence the growth of carcinoma cells 
[27,28]. Unequivocal evidence that paracrine FGF 
released from breast tumor stroma functions to promote 
tumorigenesis is lacking, but, intriguingly, we found a 
strong correlation between FGFR2 and FGF 10 mRNA 
expression levels in the cultured skin fibroblasts. We 
have not investigated the mechanism underlying this 
correlation, nor have we been able to establish that it 
also holds true for tumor-derived fibroblasts. Since 
FGF10 is known to be secreted by fibroblasts and to 
bind specifically to the FGFR2-IIIb isoform expressed on 
epithelial cells, however, our findings would fit a model 
of paracrine tumor-stroma interaction. The FGFR2-IIIb- 
FGF10 interaction plays a key role in the normal embry- 
ological and postnatal development of the mammary 
glands in mice [29,30]. Mice deficient in FgflO or Fgfr2b 
fail to develop normal mammary glands. Apparently, 
Fgfr2b signaling is crucial for the survival and prolifera- 
tion of the mammary luminal epithelial cells but does 
not affect the regenerative potential of the mammary 
epithelial progenitor cells. FgflO overexpression in the 
stromal compartment of the murine prostate results in 
epithelial cell hyperproliferation [31]. Also, in humans, 
FGF10 is thought to stimulate epithelial cell prolifera- 
tion [32]. 

Thus, it is conceivable that, in the human breast, 
slightly increased levels of FGFR2 and FGF10 result in a 
slightly increased ductal branching. If this primarily 
involved the luminal epithelial component of the breast 
tissue as it does in mice, this would explain why the 
increased breast cancer risk conferred by FGFR2 intron 
2 SNPs is mostly restricted to ER-positive tumors 
[22,33]. High FGF2 and/or FGFR2 protein levels have 
been correlated with high ER levels in breast cancer 
[34,35]. Our cohort was probably underpowered to 
demonstrate a correlation between ER status of the 
tumor and FGFR2 mRNA expression levels in skin 
fibroblasts. FGF7, also secreted by breast fibroblasts, has 
similarly been suggested to act as a paracrine growth 
factor in human breast cancer [36], but our data do not 
demonstrate a significant link between FGF7 and FGFR2 
expression. 

In pancreatic cancer, FGF10 is found in stromal cells, 
close to the tumor cells, and is thought to interact with 
FGFR2-IIIb on the tumor cells, thereby inducing cell 
migration and invasion [37]. Since the downstream 
effects of activation of FGFR2 by different FGFs in dif- 
ferent cell types are very diverse and could include 
apoptosis, cell proliferation and angiogenesis [11], it is 
difficult to predict which of these effects could mechan- 
istically explain the increased breast cancer risk 
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associated with intron 2 SNPs. Here we have explored 
some of the end points of these processes in breast 
tumors from patients with different FGFR2 mRNA levels 
in their fibroblasts. Tsunoda et al. [38] reported 
increased T-lymphocyte and macrophage infiltration 
into a newly inoculated tumor when they injected FGF2 
into this tumor. In mice, overexpression of FgflO leads 
to highly vascularized tumors in immunocompetent 
mice [39]. However, we found no correlation between 
FGFR2 mRNA levels and, respectively, stroma percen- 
tage of the tumor, microvessel density and T-cell infil- 
trate. Our sample size may have been too small to 
detect any existing differences, in particular because 
these differences can be expected to be small, given the 
mildly increased breast cancer risks conferred by FGFR2 
SNPs. 

A limitation of our study is that our observations were 
based on cultured cells from surgically removed tissues. 
Culture conditions between fibroblasts and epithelial 
cells were necessarily different and may have influenced 
some of our findings, but this enabled us to study RNA 
expression levels by cell type instead of by heteroge- 
neous tissue. It has long been known that the phenotype 
of fibroblasts can differ depending on anatomical site 
[40-43]. Indeed, we did observe absolute differences in 
FGFR2 mRNA expression levels between fibroblasts 
from skin tissue and those from normal breast tissue, 
but we demonstrated that these correlate very well (Fig- 
ure 2). 

Conclusions 

In conclusion, it is likely that the causative variants 
tagged by rs2981578 in intron 2 of the FGFR2 gene 
cause a higher breast cancer risk by influencing FGFR2 
expression levels. Here we have shown that the effect of 
intron 2 SNPs on FGFR2 expression is cell type-specific. 
Different effects were observed in skin fibroblasts and 
epithelial cells. Tissue specificity for expression quantita- 
tive trait loci is common [44] and clearly also applies to 
FGFR2. 

In addition, we observed differences in the levels of 
FGFR2 mRNA expression between fibroblasts derived 
from normal breast tissue or from tumor tissue, which 
strongly supports a holistic model to explain FGFR2- 
related breast cancer risk [45] rather than one assuming 
cell autonomous effects of FGFR2 expression modula- 
tion. Our finding that individuals in whom intrinsically 
higher levels of FGFR2 are expressed in their skin fibro- 
blasts also have higher levels of FGF10 suggests that the 
increased breast cancer risk might be due to a stronger 
paracrine effect between stromal and tumor cells invol- 
ving FGF10 signaling. The fact that the association 
between FGFR2 and FGF10 expression was not observed 
by Meyer et al. [7] in 45 normal breast tissue samples 



underscores the importance of analyzing the various 
constituting cell types separately in sufficiently large 
cohorts. Since we limited our search for such associa- 
tions to FGFs secreted by fibroblasts, and given that the 
FGFR2 and FGF10 genes are located on different chro- 
mosomes, it will be important to extend these analyses 
to genome-wide expression differences between different 
cell types in various FGFR2 genotype backgrounds. 

Additional material 



Additional file 1: Figure SI. (A) Typical fibroblast and (B) typical 
epithelial cell cultures (original magnification, xlO). 

Additional file 2: Supplementary methods This file contains 
supplementary methods regarding the quantitative real-time PCR and 
Western blot analyses. 

Additional file 3: Figure S2. Phosphorylation levels of downstream 
targets of FGFR2 after different periods of stimulation with FGF2 in one 
fibroblast sample (sample from experiment 1 with high FGFR2 mRNA 
levels). (A) Using one of the blots as an example, we show 
phosphorylated ERK1/2 examined after 0, 1, 5, 10 and 20 minutes. H: 
lanes with total protein from a fibroblast sample with high FGFR2 mRNA 
level; L: lanes from a fibroblast sample with low FGFR2 mRNA level; M: 
lane with the length marker; T: tubulin, El/2 ERK1/2. (B) Phosphorylation 
levels of FRS2a, ERK1 and ERK2 are shown as fold increases compared to 
the phosphorylation levels at 0 minutes of stimulation. 

Additional file 4: Figure S3. The relationship between FGF10 mRNA 
expression and the rs2981578 genotype in 68 skin fibroblast cultures (P = 
0.06; one-way ANOVA). The expression levels were log 2 -transformed and 
normalized to HNRPM and TBP expression. 

Additional file 5: Figure S4. Correlation of expression levels of FGFR2 
and FGF10 mRNA in 44 breast fibroblast cultures. FGFR2 and FGF10 
expression levels were normalized to HNRPM and TBP and log 2 - 
transformed. Each dot represents the expression levels of one patient 
(Spearman's p = 0.25, P = 0.11). 

Additional file 6: Figure S5. Western blot analysis measuring FGFR2 
protein levels in fibroblasts. (A) Western blot showing results for seven 
different fibroblast cultures. H: lanes with total protein from fibroblasts 
with high FGFR2 mRNA levels; L: lanes with total protein from fibroblasts 
with low FGFR2 mRNA levels; empty: empty lane; M: lane with length 
marker; recomb: lane with pure recombinant FGFR2. (B) FGFR2 mRNA 
levels and FGFR2 protein levels in the seven fibroblast samples. The 
FGFR2 mRNA results were normalized to HNRPM and FBP and log 2 - 
transformed. The FGFR2 protein results were normalized to a-tubulin. 
Each dot represents the results for one fibroblast sample. 
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